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Abstract. Single-chain magnets are molecular spin chains displaying slow relax- 
ation of the magnetisation on a macroscopic time scale. To this similarity with 
single-molecule magnets they own their name. In this chapter the distinctive fea- 
tures of single-chain magnets as opposed to their precursors will be pinpointed. In 
particular, we will show how their behaviour is dictated by the physics of thermally- 
excited domain walls. The basic concepts needed to understand and model single- 
chain magnets will also be reviewed. 



1.1 Introduction 



The observation of magnetic hysteresis of molecular origin in Single-Molecule 
Magnets (SMMs) is considered one of the most relevant achievements in nano- 
magnetism [1,2]. Fundamental aspects related to quantum tunnelling of the 
magnetisation have been thoroughly discussed in the previous chapters. On a 
more practical perspective, that observation rendered the molecular approach 
one of the possible routes to realizing bistable nano-objects, suitable for mag- 
netic storage or quantum-computing applications. In spite of many efforts, 
the highest blocking temperature attained by SMMs remains, still nowadays, 
in the liquid-helium temperature range. The idea that one-dimensional (ID) 
structures of coupled paramagnetic ions might afford higher blocking temper- 
atures started developing at the end of the nineties and the first examples of 
slowly relaxing ID systems were reported at the beginning of the new cen- 
tury [3,4]. The resulting molecular systems have been dabbed Single-Chain 
Magnets (SCMs) in order to evidence analogies with their precursors, SMMs, 
while remarking - at the same time - the ID character. In some cases, SMMs 
themselves have been employed as building blocks for such ID magnetic lat- 
tices [5] . With the aim of increasing the blocking temperature as much as pos- 
sible, different synthesis strategies have been followed to obtain some type of 
magnetic anisotropy at the level of building blocks or of the coupling among 
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thorn. In the present chapter we will be dealing with uniaxial anisotropics 
only, though this requirement is not strict for the observation of SCM be- 
haviour [7-9]. 

A distinctive feature of ID magnetic systems is the development of short- 
ranged correlations upon cooling. This makes them substantially different 
from both paramagnets and bulk magnets. Should one establish an analogy 
between classical magnetic ordering and phases of matter, paramagnets would 
be identified with perfect gases while bulk magnets with solids. Pushing this 
naive analogy further, spin chains would be associated with liquids, specifically 
in the temperature range in which short-range correlations extend over several 
lattice units. The degree of spatial correlation is quantified by the correlation 
length. In molecular chains consisting of magnetic building bocks with uni- 
axial anisotropy, the correlation length typically diverges exponentially with 
decreasing temperature. From a snapshot taken at finite temperature, any 
chain would appear as a collection of randomly oriented magnetic domains 3 
separated by domain walls (DWs). The average size of those domains is of 
the order of the correlation length. This pictorial, but essentially correct, sce- 
nario is consistent with thermally-driven diffusion of DWs. In this sense, the 
response of a SCM to a tiny a.c. field is expected to be determined by the time 
needed to adjust the size of domains to the external stimulus. A robust scaling 
argument associates the characteristic time of this readjustment with the time 
elapsed while a DW diffuses over a distance proportional to the correlation 
length. Within this idealized picture, the relaxation time of the magnetisation 
is expected to scale with temperature like the square of the correlation length. 

The qualitative description given above applies to the ideal case of infi- 
nite chains and small applied fields. The first hypothesis is practically never 
fulfilled in real systems. In fact, the number of magnetic centres interacting 
consecutively is typically limited to 10 2 -10 4 by naturally occurring defects, 
non-magnetic impurities or lattice dislocations [10-12]. A SCM may thus be- 
have as if it extended indefinitely only when the correlation length is much 
smaller than the average distance between successive defects. Upon lowering 
the temperature, a crossover is expected at which the correlation length be- 
comes of the order of the average distance among defects. Below this crossover 
temperature, spins enclosed between two successive defects are parallel with 
each other and no DW is present at equilibrium. In this finite-size regime re- 
laxation is somewhat equivalent to magnetisation reversal in nanoparticles or 
nanowires, which may occur via Neel-Brown uniform rotation or by droplet- 
nucleation mechanism [13]. 

All the mentioned mechanisms represent potential channels for relaxation 
in SCMs. Which one is favoured depends on the experimental conditions: tem- 
perature, applied field and amount of defects in the sample. The random-walk 
argument which relates the correlation length to the relaxation time holds in 



3 These soft, fluctuating domains should not be confused with Weiss domains en- 
countered in magnetically ordered phases. 
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the linear-response regime, i.e., when such tiny fields are applied to induce 
just slight deviations from thermodynamic equilibrium. On the contrary, re- 
laxation from a saturated configuration typically entails far-from-equilibrium 
dynamics. In this type of experiments nucleation of soliton-antisoliton pairs 
or of a single DW adjacent to a defect possibly initiates the relaxation pro- 
cess. Neel-Brown uniform rotation practically represent an alternative channel 
for relaxation only for very short segments of chain, encountered in samples 
in which finite-size effects have been enhanced by doping with non-magnetic 
impurities [10, 11]. 

It should not be forgotten that molecular spin chains are packed in three- 
dimensional crystals. Though several synthesis strategies may be followed to 
minimize interactions among chains, at least the dipolar interaction cannot be 
suppressed completely. Therefore, below some temperature, a 3D magnetically 
ordered phase is expected to appear. Whether such a phase is observed or not 
in a specific compound depends on how long the relaxation time is at the 
transition temperature [14]. Generally, when the time needed for the system 
to equilibrate is much longer than experimental time scales, the distinctive 
features of the underlying equilibrium phase, possibly ordered, cannot be ev- 
idenced. For weakly interacting spin chains, the expected transition tempera- 
ture is much higher than interchain interaction in Kelvin units (fcTc 3> J' with 
the forthcoming notation). In fact, the 3D-ordering process is "assisted" by 
the development of strong short-range correlations inside each chain [15-17]. 
However, in realistic samples, defects prevent the intrachain correlation length 
from diverging indefinitely, which eventually lowers the transition tempera- 
ture to the ordered phase. In several SCMs slow dynamics was observed down 
to few Kelvins, before 3D ordering took place, right because of the presence 
of defects and non-magnetic impurities. 

Both SMMs and SCMs are characterized by slow dynamics of molecular 
origin, acting at macroscopic time scales and in the absence of 3D magnetic 
ordering. Even if impurities play a crucial role in SCMs, they usually do 
not bring enough disorder to give rise to spin-glass behaviour. Consistently, 
slow dynamics is typically characterized by a single time scale which does 
not display a super- Arrhenius behaviour at any temperature [18]. Besides 
preventing the onset of 3D magnetic ordering, the increase of relaxation time 
with cooling usually leads to complete blocking before genuine quantum effects 
become evident [19]. 

From what written till now, it should be clear that many effects may 
interplay in determining the magnetic behaviour of spin chains. We will focus 
on those systems in which slow dynamics can be ascribed to each single chain 
and does not originate from cooperative 3D interactions. 

The goal of this chapter is that of highlighting the properties of SCMs 
with a critical view to what has been done and what still deserves further 
investigation. We will not try to cover in detail all the representative litera- 
ture, for which the reader is addressed to excellent reviews [6-9]. Although 
SCMs have been widely investigated, the interest of the physics community 
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has not been comparable to that shown for SMMs. This is partially due to 
the fact that the novelty of SCMs compared to traditional ID spin systems 
hardly emerged. Here we attempt to provide and efficient overview of the es- 
sential, novel physics of SCMs and hope that the final comment will be more 
benevolent. The chapter is organised as follows: Sect. 1.2 will cover the basic 
aspects of classical spin chains; the chemical frame will be discussed within 
a bottom-up or building-block approach in Sect. 1.3; in Sect. 1.4 the spin 
Hamiltonians typically used to rationalise the physical properties of SCMs 
will be introduced; relevant extensions of the Glauber model developed in the 
SCM context without and with defects will be treated in Sect. 1.5 and 1.6, re- 
spectively; phenomenological arguments not contained in the Glauber model 
but relevant for understanding SCMs will be discussed in Sect. 1.7; a short 
section on perspectives will conclude the chapter. 

1.2 Thermal equilibrium and slow dynamics in ideal 
SCMs 

In this section the peculiarities of classical spin chains with uniaxial anisotropy 
that directly affect the physics of SCMs will be recalled. Indeed, the distinctive 
feature of SCMs is that of approaching thermodynamic equilibrium slowly. 
By slowly we mean that relaxation time becomes longer than milliseconds 
at temperatures of the order of 10 K or lower. The reference equilibrium 
state to be reached is also relevant. As already mentioned, as long as 3D 
interactions are negligible, no magnetisation is expected in zero applied field at 
thermodynamic equilibrium. Long-range magnetic order may be destroyed by 
thermally-excited spin waves or DWs either. The first ones are effective in the 
absence of anisotropy, according to the Mermin-Wagner theorem [20,21]. The 
fact that disordering is, instead, driven by DWs in the presence of anisotropy 
can be easily understood recalling an argument presented in the Landau- 
Lifshitz series [22] . Let us consider a group of N spins that preferentially point 
along the same direction, say up or down. For the moment we assume the axes 
of easy anisotropy to be collinear, as represented schematically in Fig. 1.1. We 
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Fig. 1.1. Sketch representing the configurations whose free-energy difference is eval- 
uated in the text: a ferromagnetic ground state with all the spins parallel to each 
other (right) and a configuration consisting of two domains with opposite spin align- 
ment (left). 
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evaluate the variation of the free energy associated with the creation of a DW 
starting from a configuration with all the spins parallel to each other. Creating 
a DW increases the energy by a factor E 2 — E\ = £dw On the other hand, 
such a DW may occupy N different positions in the spin chain, so that the 
relative entropy increase scales as 52— Si — k ln(iV). The free-energy difference 
between the two configurations sketched in Fig. 1.1 is roughly AF ~ £dw — 
kTln(N). When the thermodynamic limit N — > oo is taken, one immediately 
realizes that it is always convenient to split the system into groups of parallel 
spins. As a consequence, long-range magnetic order is destroyed at any finite 
temperature. In principle, in an infinite chain, the same mechanism may allow 
creating an indefinite number of DWs. However, the average distance among 
them does depend on temperature and it is inversely proportional to the 
correlation length [23]. It is worth remarking that in the text-book argument 
given above the following assumptions have been implicitly made: 

1. that DWs extended only over one lattice unit 

2. spin-wave excitations were not considered 

3. the thermodynamic limit was taken. 

Whether the first hypothesis is fulfilled or not depends on the relative strength 
of exchange interaction and magnctic-anisotropy energy. This can be discussed 
more concretely by considering the classical Hciscnbcrg model with uniaxial 
anisotropy: 

N 

n H = -J2[ Js ^- s i+1 + D(s?f] , (i.i) 

where Si are classical spins, J and D the exchange and the magnctic- 
anisotropy energy, respectively; \Si\ = 1 will be assumed henceforth. Though 
it does not entail the complexity of many real systems, Hamiltonian (1.1) 
is a useful reference to discuss the physics of SCMs. To the aim of distin- 
guishing between two simple types of DWs, we fix D > and J > 0. With 
Hamiltonian (1.1), DWs can be larger than one lattice spacing. In fact, the 
actual DW profile results from the competition between the exchange energy 
(which is minimized by broadening the wall) and the anisotropy energy (which 
favours a sharp wall) . Domain walls whose structure develop over more lattice 
units will be referred as broad; these are opposed to sharp DWs in which the 
local magnetisation changes abruptly its sign, within one lattice distance. The 
energy associated with a broad DW is £dw = 2v / 2~DJ [24], namely the energy 
needed to create one soliton "particle" in the spin chain [25]. For sharp DWs, 
one obtains £dw = 2J, as per the Ising model. The crossover from sharp- 
to broad- wall occurs at D/J = 2/3 [26-28]. The analytic formula for broad- 
DW energy, £ d w = 2y/2DJ, was obtained in the continuum formalism and 
gets less and less accurate as the threshold ratio is approached from below, 
D/J -> (2/3)-. 

If the Landau's argument is rephrased for DW excitations of finite thick- 
ness w = yj J/2D, the counting of equivalent configurations with the same 
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energy needs to be modified and - in turn - the entropy contribution 
S2 — Si = kln(N/w). In this case, splitting the uniform configuration into 
domains becomes convenient as soon as the number of spins exceeds the 
product w e £dw / fcT . The latter threshold gives an estimate of the average num- 
ber of consecutive spins that can be found aligned at a given temperature. 
To the leading order, the correlation length scales in the same way at low 
temperature: £ ~ i»e £d, ' tT . The energy £dw represents the natural "unit" 
which controls the divergence of the correlation length. Thus, in classical spin 
chains with uniaxial anisotropy the characteristic exponential divergence of £ 
is closely related to the fact that ferromagnetism is destroyed by thermally 
excited DWs. 
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Fig. 1.2. Log-linear plot of £ in lattice units computed with the transfer-matrix 
technique as function of the ratio £d v /kT for different values of D/J. For D/J= 
0.1 (red crosses), 0.3 (green crosses), 0.5 (blue stars) DWs are broad and £d w has 
been computed numerically on a discrete lattice. For D/J— 5 (open squares) DWs are 
sharp and £dw=2 J has been used. The two solid lines give the "reference" behaviour 
£~e /1 £' /feT which is indeed followed when DWs are sharp (D/J=5) but not when 
DWs broaden. Inset: A=£/w is plotted as a function of S^/kT for the values of 
D/J consistent with broad DWs. The universality of A is highlighted by the data 
collapsing. Solid lines evidence the decrease of A% with increasing temperature [30] . 

In contrast to the Ising model [29], the classical spin Hamiltonian (1.1) 
can also host spin-wave excitations, besides DWs. Due to the interaction be- 
tween spin waves and broad DWs an additional temperature-dependent factor 
appears in front of the exponential in the low-temperature expansion of the 
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correlation length [31]. Moreover, spin waves renormalise the DW energy at 
intermediate temperatures. The net result of the complicated interplay be- 
tween thermalised spin waves and DWs is that the energy barrier controlling 
the divergence of £ (usually called A$ in SCM literature [6-9]) is generally 
smaller than £dw and takes different values depending on the temperature 
range in which it is measured [30]. A similar effect was reported for the ac- 
tivation energy of 2n sine-Gordon solitons in Mn 2+ -radical spin chains [32]. 
Fig. 1.2 highlights how A^ is constant and equal to £dw for sharp DWs, while 
it varies significantly for broad DWs. However, the correlation length in units 
of w keeps depending only the ratio Sdw/kT, i.e., £/w = A (£& w /kT). The 
inset shows that the curves corresponding to broad DWs indeed collapse onto 
each other when the ratio £/w is plotted as a function of £& w /kT. 

As mentioned in the introduction, in realistic spin chains the divergence 
of the correlation length is always hindered by the presence of defects and 
non-magnetic impurities. This implies that results derived taking the thermo- 
dynamic limit, N — > oo, do not hold down to indefinitely low temperatures. If 
we assume - for the time being - an idealized scenario in which such defects 
do not occur, a certain number of DWs shall be present at any finite tem- 
perature. A simple random- walk argument then relates the relaxation time to 
the correlation length: within a time r a DW performs a random walk over a 
distance proportional to £ [33]. In other words, the relation 

£ 2 * 2D s t (1.2) 

can be assumed, with D s being the diffusion coefficient. The latter generally 
increases with increasing temperature. Moreover, it is expected to depend on 
temperature differently for sharp or broad diffusing DWs. When presenting 
the Glauber model we will see that D s can also be interpreted as the attempt 
frequency to flip a spin adjacent to a sharp DW. 

Summarizing, the presence of uniaxial anisotropy produces an exponential 
divergence of the correlation length with decreasing temperature. As the re- 
laxation time is related to £ by a random-walk argument, it is also expected 
to diverge likewise, so that 

^e A ^ kT r~e^/ fcT . (1.3) 

In ideal ID magnetic systems [25,34] the correlation length is proportional to 
the product of temperature by static susceptibility (measured in zero field): 

XcqT~£. (1.4) 

The relaxation time can, instead, be obtained from dynamic susceptibility 
measurements as follows 
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where oj is the frequency of the oscillating applied field and Xoq is the static 
susceptibility 4 . Both real and imaginary part of x{ UJ ^ T) display a maximum 
for lot = 1. The basic experimental characterization of SCMs essentially re- 
duces to determining the temperature dependence of £ and r, which is - in 
principle - possible thanks to equations (1.4) and (1.5). 

Even within the idealized scenario presented in this section, the way in 
which and A T defined in (1.3) relate to the Hamiltonian parameters J and 
D depends on the DW thickness, w, and on the temperature range in which 
such energy barriers are measured. Besides this, model Hamiltonians of real 
SCMs may differ significantly from (1.1). In the next section we will recall 
some synthesis strategies that have been followed to produce different SCMs. 
The features of the employed building blocks and the type of coupling among 
them eventually decide which model is more appropriate to describe a specific 
SCM. 

1.3 Tailoring SCMs by building-block approach 

The initial interest on Molecular Magnets stemmed from the attempt to de- 
sign molecular systems displaying long-range magnetic order. However, after 
more than 30 years of attempts there are only two room-temperature molec- 
ular magnets and matters are no better for liquid- nitrogen temperatures [36] . 
To have long-range order it is necessary to build 2D or 3D structures of centres 
magnetically coupled. This is difficult with molecular bricks since the number 
of coordination sites which are available to propagate the exchange coupling 
in different directions is small due to the presence of capping ligands. Such 
bricks are then more suitable to produce low-dimensional systems, like clus- 
ters of metal ions (zero dimensional) [2] or spin chains. These systems do not 
display long-range order but still show a variety of interesting phenomena, 
including SCM behaviour. It is pedagogically useful to imagine that synthe- 
sizing a SCM is like assembling bricks with a magnetic functionality and a 
structural functionality. Usually, the latter is provided by organic molecules 
and the former by metal ions. Building blocks need to be chosen and arranged 
in a structure which maximizes the intrachain and minimizes the interchain 
interactions. Bricks with magnetic functionality must be coupled ferro- or 
ferri-magnetically and control of the magnetic anisotropy must be achieved. 
Chemists are not yet able to have that detailed control but progress is fast 
and serendipity always helps. 

Some centres of the building blocks shall be magnetically active, which 
implies the presence of unpaired electrons that are formally assigned to mag- 
netic orbitals, either p, d, or f. In organic radicals the unpaired electrons 
normally belong to p orbitals: these are external orbitals which strongly in- 
teract with the environment. For this reason such electrons hardly remain 

4 The more general Cole-Cole equation is needed when relaxation is not charac- 
terised by a single r or to account for adiabatic contribution to x [35] . 
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unpaired but rather tend to couple with electrons of other molecules in cova- 
lent bonds, which eventually explains why few stable organic radicals exist. 
In the following, we will mostly refer to nitronyl nitroxide radicals (NITR), 
whose structure is shown in Fig. 1.3a. The unpaired electron is delocalised on 
the group O-N-C-N-0 and, from the magnetic point of view, basically behaves 
as a free electron. Its magnetic moment is essentially spin determined, with 
little orbital contribution due to small spin-orbit coupling. This implies low 
magnetic anisotropy which is the final blow for purely organic SCMs. 

NITR trans cis 



group / 

\ i 

c c 



(a) (b) (c) 

Fig. 1.3. (a) schematic structure of the NITR radical: the unpaired electron is 
delocalised over the O-N-C-N-O fragment (coloured), which is magnetically active, 
(b) and (c) show two possible realisations of direct exchange coupling between the 
electron of each NITR radical and a metal ion (intrachain interaction): each M can 
be bond to two NITR groups through oxygens occupying either trans (b) or cis (c) 
positions in the coordination octahedron (see Fig. 1.4). 

NITR radicals have the right geometry for bridging two metal ions through 
their equivalent oxygen atoms (extremes of the O-N-C-N-0 fragment in 
Fig. 1.3a). The above considerations suggest that NITR radicals are not ap- 
propriate for being used alone, but they become excellent bricks for SCMs 
when coordinated to metal ions [6-9,37]. In fact, the interaction of the p or- 
bitals with the d (or f) orbitals can be strong, of direct type, both ferro and 
antiferromagnetic in nature. 

Transition-metal ions provide good magnetic bricks. As anticipated in the 
introduction, we will limit ourselves to consider SCMs possessing uniaxial 
anisotropy at the brick level. In molecular systems, magnetic anisotropy is 
closely related to the fact that the surrounding of metal ions is not spherically 
symmetric. Figure 1.4a shows a generic metal atom (M) in an octahedral en- 
vironment of ligands. Oxygen atoms occupy the vertices of the octahedron. In 
the group M(hfac)2, for instance, two oxygens of each hexafluoroacetylaceto- 
nate (hfac) ligand coordinate to M, thus occupying two neighbouring vertices 
of the octahedron per hfac molecule. The two remaining, empty coordination 
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sites can be in either trans or cis position (Fig. 1.4 b and c, respectively) 
and may host oxygens of other ligands that can be used to connect differ- 
ent M(hfac)2 moieties. The choice of NITR to bridge those moieties creates 
a strong, direct exchange coupling between M and the electron delocalised 
on each O-N-C-N-0 group (intrachain interaction, J). Consistently with the 
two possible coordination configurations of M(hfac)2 sketched in Fig. 1.4, the 
segments connecting the metal ions with NITR oxygens may form an angle of 
180° (trans) or 90° (cis). Such segments specify the direction along which the 
intrachain interaction propagates. The bulky hfac groups prevent efficient in- 
terchain exchange coupling. The residual interchain interaction J' has mainly 
dipolar origin and it is, typically, from 3 to 6 orders of magnitude smaller 
than the intrachain interaction. For this reason, M(hfac)2 moieties are per- 
fectly suited for realizing isolated spin chains (| J' / J\ < 10~ 3 ). The interaction 
between successive magnetic bricks can either be ferro- or ferrimagnetic and 
give rise to straight or zig-zag structures. Besides, magnetic bricks are often 
characterized by low symmetry with the metal ions occupying general posi- 
tions in the unit cells, which does not impose limitations to the orientation 
of anisotropy axes. Therefore, in practice, full collinearity among anisotropy 
axes is more an exception rather than the rule [6-9]. 




(a) (b) (c) 



Fig. 1.4. Sketch of a metal ion (purple spheres) in an octahedral environment of 
oxygen- donating ligands (red spheres representing oxygens), (a) metal-oxide coor- 
dination in an extended solid, (b) and (c) M(hfac)2 moiety with two empty coor- 
dination sites in trans position (b) and in cis position (c); the CF3 groups of hfac 
ligands are not shown for clarity sake. Green dashed lines indicate the directions 
along which the intrachain exchange coupling mediated by a different ligand (e.g., 
NITR radical) may propagate. 

Since SCM behaviour requires some magnetic anisotropy, the orbital mo- 
mentum must not be completely quenched. The surviving component may 
be associated with single-ion anisotropy or with pair-spin interaction. In the 
former case, the residual orbital contribution can show up in a g tensor differ- 
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ent from the free-electron one and/or in the zero-field splitting. With a large 
periodic table it is amazing that only cobalt and manganese, with some iron 
and nickel have been used. Mn 3+ is an example of anisotropy determined by 
zero-field splitting; while in Co 2+ the anisotropy is associated with the g ten- 
sor [38] . The crystal- field theory is the simplest way to describe the ground and 
low-lying levels of a transition-metal ion. The Hamiltonian can be expressed 
as a sum of terms: 

H = Ho+n cc + H CF +H LS , (1.6) 

where "Ho is the origin of the electron configuration (3d)™, H cc is the electron- 
electron repulsion, Hqf is the crystal-field term and Hls is the spin-orbit cou- 
pling. For 3d ions "Hcc and "Hcf are comparable and larger than spin-orbit cou- 
pling. It is customary to neglect in first approximation the spin-orbit coupling 
which is introduced later as a perturbation. Mn 3+ has a (3d) 4 valence-electron 
configuration which in octahedral symmetry yields a 5 E g ground state. This 
is unstable to phonon coupling (Jahn- Teller theorem), which lowers the sym- 
metry to D4h, namely to a tetragonally elongated coordination. The ground 
state 5 Ai g , in zero-order approximation, is five-fold degenerate (no orbital 
degeneracy and 5=2 spin multiplct). The spin-orbit perturbation yields no 
contribution in the first order, but to the second order it admixes excited 
states with the ground state. This removes the degeneracy of the spin multi- 
plet and produces anisotropic components in the g tensor. Its effect is usually 
summarized introducing an effective single-ion spin Hamiltonian of the form: 

H = -DS 2 Z - fi B B g S , (1.7) 

where g is a symmetric tensor. The first term is responsible for the zero-field 
splitting of the 25+1 levels. It is often referred to as crystal-field term, even 
though this is misleading because it is not the crystal field which splits the 
levels but rather the spin-orbit coupling. The spin-Hamiltonian parameters are 
determined by the spin-orbit coupling constant A and by the degree of mixing 
between the 5 Ai g ground state and the excited states induced by "Hls [39]. 
The lowering of the symmetry produces axially symmetric g and D tensors 5 
whose components are related, to the leading order, through the following 
formula: 

D = D Z - D x . y = X {g Xiy - 9z ) . (1.8) 

1.4 Realistic spin Hamiltonians for Single-Chain 
Magnets 

So far we have neglected the coupling among spin pairs, which can be written 

as 

5 Without loss of generality g and D can be assumed symmetric. Consequently, 
they are diagonal on a proper reference frame with eigenvalues g x , g y , g z and D x , 
D y , D z . The same notation will be used for the J tensor. 
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"Hcxch = —S p JSk , (1.9) 

where S p and Sk are effective spin operators of any two interacting magnetic 
bricks and J is a generic 3-by-3 matrix. Limiting oursclf to intrachain spin-spin 
coupling, we can neglect the contribution due to dipolar interaction which is 
typically much smaller than the exchange one. When pair-spin interaction in- 
volves transition metals whose ground state is not orbitally degenerate, the 
isotropic contribution to the J tensor dominates. As mentioned before, second- 
order perturbation theory prescribes that the ground-state wave functions be 
modified because of the admixing with excited states mediated by spin-orbit 
coupling. When the corrected wave functions of the bricks p and k are em- 
ployed to compute the exchange integral, the anisotropic and antisymmetric 
contributions emerge. The former is proportional to (Ag/g) 2 , while the latter 
is proportional to Ag/g. This ratio is usually much smaller than one, thus the 
antisymmetric term - if allowed by symmetry - is expected to dominate with 
respect to the anisotropic exchange. For our purposes, it will be enough to 
know that anisotropic contributions to J can be neglected when the g tensor 
is fairly isotropic, as for Mn 2+ , Mn 3+ , high-spin Fe 3+ , etc. 

When magnetic bricks comprise transition metals with orbitally degener- 
ate ground state, predicting the properties of the g, D and J tensors on simple 
footing becomes extremely complicated [40] . One possibility is that of consid- 
ering just a symmetric exchange tensor, which is then diagonal on a proper 
basis with principal values J Xl J y and J z . If compatible with symmetry, an 
antisymmetric, Dzyaloshinskii-Moriya term may be added independently. 

In ID magnetic systems realized by coupling radicals with neighbouring 
transition-metal ions, anisotropic terms in J may originate only from the metal 
atoms. The first successful examples consisted in ferrimagnetic chains of gen- 
eral formula Mn(hfac)2NITR [17]. The radical is isotropic and so is Mn 2+ , 
therefore J is expected to be proportional to the identity. Indeed, these sys- 
tems represented text-book examples of ID Hcisenberg ferrimagnets described 
by the Hamiltonian 

N/2 

'HMn-rad = ~ J ^ S 2 p ■ S 2p+ \ , (1-10) 
p=l 

where S 2p stand for Mn 2+ spin operators (lying at even sites 2p with 
S 2p = 5/2), while s 2p+ i are the radical spin-one-half operators. J is nega- 
tive and tends to orient the nearest-neighbouring spins antiparallel to each 
other. The temperature dependence of the static susceptibility was fitted us- 
ing the Seiden model [41] with \J\ in the range 300-475 K depending on the 
substitucnt R on the radical 6 [42]. In the Seiden model the Mn spins are 
replaced by classical vectors, which - in the absence of field and single-ion 

6 Henceforth, energies will be expressed in Kelvin units to make it easier to com- 
pare them with thermal energy. The conversion factor to SI coincides with the 
Boltzmann constant k: 1 K = 1.3806503 x 10 -23 J. 
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anisotropy - makes the model analytically solvable. Due to the large value 
of the coupling between Mn 2+ and NITR radicals, strong pair-spin correla- 
tions develop, which is highlighted by a divergence of the correlation length 
at low temperature (proportional to \ J\/T and not exponential like in spin 
chains with uniaxial anisotropy [42]). In the presence of such strong intra- 
chain correlations even a tiny interchain interaction J' may induce 3D or- 
dering [15]. In Mn(hfac) 2 (NITiPr) this happens at T c =7.6 K [17]. ESR and 
NMR studies provided evidence of spin-diffusion effect allowing for an esti- 
mate of the ratio between inter- and intrachain exchange interaction of the 
order |J'/J|=2xlO- 6 [16]. 

This example confirms that combining transition metals with organic radicals 
is a powerful strategy for designing ideal ID systems. An additional ingre- 
dient is needed to realize a SCM: magnetic anisotropy. This may easily be 
introduced by replacing Mn 2+ with Mn 3+ . Recently, the observation of slow 
relaxation consistent with SCM features was reported for ferrimagnetic spin 
chains consisting of Mn 3+ and TCNQ or TCNE 7 organic radicals [9,43]. The 
relatively large multiplicity of Mn 3+ spins, £=2, allows justifying their re- 
placement by classical vectors. Thus, the Seiden model is still a good starting 
point for describing the magnetic properties of these systems, provided that 
single-ion-anisotropy terms are added. Even if the modelling aspects are well- 
defined, the rationalization of Mn 3+ -radical SCMs is complicated by the fact 
that J 3> D, meaning that the relevant excitations are broad DWs. 

The extreme anisotropic g tensor obtained for Co 2+ in a tctragonally com- 
pressed symmetry suggests that its coupling with NITR be, to leading order, 
of the Ising type. This idea led to the synthesis of the first compound showing 
SCM behaviour: Co(hfac)2(NITPhOMe) [3]. Experimental results pertaining 
slow dynamics have shown a substantial agreement with the kinetic version 
of the Ising model developed by Glauber [44] . Unfortunately, up to date, the 
static properties have not been successfully modelled yet. The first reason is 
that above 40 K treating Co 2+ as an effective S =1/2 is not legitimate (the 
energy separation between the ground-state Kramers doublet and the excited 
multiplets is about 100 K). A second reason relates to the helical structure of 
this compound, because of which the elementary magnetic cell contains 3 Co 2+ 
and 3 radical spins. Apart from the question of reproducing its static proper- 
ties, it is instructive to give a closer look at the Hamiltonian of this system to 
show how non-collinearity can be modelled in general. For temperatures lower 
than 40 K, a reasonable Hamiltonian for the Co(hfac)2(NITPhOMe) chain is 
given by 

N/N r N r /2 

"Hco-rad = ~ ^S p ,2r h r Sp,2r+l + A*B-Bg2r S Pi 2r + gB ■ S p , 2r +1 , 

p— 1 r— 1 

(1.11) 

7 Acronyms stand for tetracyanoquinodimethane (TCNQ) and tetracyanoethylene 
(TCNE). 
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where both S Pt2r and s p ^ r +i are spin one-half operators associated with Co 2+ 
ions and radicals, respectively, p represents the magnetic cell index while r 
spans the inequivalent Co 2+ atoms inside each cell. For the specific case, r 
takes N r =3 different values which correspond to different orientations of the 
principal axes along which the J and g tensors are diagonal. If spin projections 
are expressed in the crystal frame, the tensors appearing in Hamiltonian (1.11) 
are built applying a standard 0(3) rotation to the diagonal tensors [38,47,48]. 
Formally, r in J 2r and g 2r labels different sets of rotation angles. The Lande 
factor of the radical is isotropic and thus independent of r. 
When spins S > 1/2 are considered, a magnetic brick may possess some single- 
ion anisotropy, which implies that also the D tensor needs to be rotated in 
non-collinear systems. 

The thermodynamic properties of classical spin chains with nearest-neighbour 
interactions can be efficiently computed by means of the transfer-matrix 
method. Letting the general Hamiltonian be H = —kT^2 p V(S p , S p+ i), the 
partition function Z is obtained integrating over all the possible directions 
along which each unitary vector S p may point: 

Z = Jdf2 1 Jdf2 2 ...J c v ^' s ^ c v ^' s ^ . . . e y ( s - Sl ) d!2 N . (1.12) 

Defining the transfer kernel as JC(S P , S p+ i) — e v ( Sp ' Sp+1 ^ and assuming pe- 
riodic boundary conditions, the partition function Z can be recasted into the 
trace of the iV-th power of JC(S P , S p+ \): 

Z = J dQ l J dfl 2 ... J K(S U S 2 )K,(S 2 , S 3 ) ... K,(S N , S{) dn N = Tr{/C N } . 

(1.13) 

When the transfer kernel is expressed on a basis of eigenf unctions, the par- 
tition function reduces to a sum of eigenvalues Z — J2 m ^Ii where ip m (S p ) 
and A m are solutions of the following eigenvalue problem: 



/ 



JC{S P , S p+1 )ip m (S p+1 )dQ p+1 = X m tp m (S p ) . (1.14) 



For kernels that can be written in a symmetric form with respect to the ex- 
change S p <-» S p+ i the spectral theorem warrants that eigenvalues are real 8 . 
They are also positive, because the transfer kernel is a positive function of S p 
and S p+ i, and upper bounded so that they can be ordered from the largest 
to the smallest one: Ao > Ai > A 2 > ... In the thermodynamic limit the 
asymptotic behaviour of the partition function (1.13) is dominated by the 
largest eigenvalue A , Z ~ A^, meaning that the free energy per spin is given 
by / = — fcTlogAo- Macroscopic experimental observables are obtained as 



In the general, non-symmetric case left and right eigenfunctions of JC(S P , S p+ i) 
have to be considered, but the basic ideas of the transfer-matrix method remain 
the same. 
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derivatives of /, but this method allows computing microscopic averages as 
well. Apart from some fortunate cases [21], equation (1.14) needs to be solved 
numerically by sampling the unitary sphere with a finite number of special 
points. This number can be increased dynamically untill the desired precision 
is reached. Even though it may not be transparent from our description, a 
new eigenvalue problem ought to be solved for any computed temperature or 
applied field. Referring the reader to the existing literature for implementation 
details [45-47] , we remark that the transfer-matrix method allows computing 
the magnetic properties of classical spin chains more efficiently than, e.g., stan- 
dard Metropolis MonteCarlo. This makes it possible to fit spin Hamiltonian 
parameters directly to experimental data sets. The major drawback is that 
the number of spin variables that appear in the kernel scales like the range 
of interaction (2 for nearest-neighbour, 4 for second nearest-neighbour inter- 
action, etc.), which finally affects the complexity of the eigenvalue problem 
in (1.14). 

The transfer-matrix method can easily be extended to models in which 
classical and quantum spins alternate, like in the Seiden model [41]. Noting 
that the quantum-spin operators are not directly coupled with each other, 
one can integrate out their degrees of freedom independently In fact, a 
generic quantum spin located at site 2p + 1 experiences an effective "field" 
kTh,2 P .2p+2 = J {^2p + S2P+2) + ^BgB. The corresponding energy levels are 
±feT|/i2p,2p+2|j which depend parametrically on the orientation of the two 
classical spins, S^p and 5^+2- After tracing over the quantum degrees of 
freedom, one is left with the kernel 

r^c c \ wim h ( mB g S 2 p \ ( D{S Z 2 P ? \ 
lC{S 2p , S 2 p+2) = 2cosh(|ft, 2P: 2p+2|)exp I fcy I exp I — I 

(1.15) 

where the single-ion anisotropy and Zeeman term acting on the classical spins 
have been added. The kernel (1.15) may be used to compute, e.g., the equi- 
librium suceptibility of Mn 3+ -radical chains [43]. To the aim of sketching how 
to proceed for modelling non-collinearity, let us substitute the spins S p ^ r m 
Hamiltonian (1.11) by classical vectors 9 . Even after integrating out the radical 
degrees of freedom, the are still 3 non- equivalent classical spins in each mag- 
netic unit cell, resulting in 3 different kernels if B is applied along a generic 
direction: /C(S P;2 , S p ^), JC(S p a, S Pi q) and K.(S P fi, S p+ i^)- The role of the 
kernel (1.15) is played by 

fc(Sp,2, Sp+1,2) = J dQi J /C(S Pj2 , S Pt i) /C(S P> 4, S P fi)1C{Spfi, S p +i j2 ) dQ§ , 

(1.16) 



Even though this is not justified for the specific case of Co 2+ , the classical ap- 
proximation allows us to discuss the general formalism. 
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Fig. 1.5. Sketch of the periodicity associated with Hamiltonian (1.11): small, black 
arrows represent radical spins s Pl 2r+i, while the large coloured arrows represent 
the metal-ion ones S Pt 2r- The kernel (1.16) depends only on S p ,2 and S p +i,2 (blue 
arrows) because an integration over the internal dgrees of freedom S Pt 4 (red arrows) 
and S p ,6 (green arrows) has been performed. 

obtained by tracing over the degrees of freedom internal to the considered cell, 
dl?4 and di?g 10 . The way in which the kernel is built is sketched pictorially 
in Fig. 1.5. In the thermodynamics limit the partition function is given by 
Z ~ Aq Nr j where the number of spins have been replaced by the number of 
unit cells N/N r . 

Due to non-collincarity, a strong anisotropy in the D, g or J tensors may not 
necessarily be evident at the macroscopic level [48]. More concretely, having 
similar saturation values for the magnetisation along different crystallographic 
directions may still be compatible with a strong uniaxial character at the level 
of individual bricks. Non-collincarity is also consistent with an inversion of the 
directions of easy and hard magnetisation by increasing temperature [47] or 
with the vanishing of the correlation length for some specific applied fields [49] . 
In passing, we note that finite-size effects can be taken into account in the 
general transfer-matrix framework [45] as well as interchain interactions if 
treated at the mean-field level [15,50]. 

In the cases in which one of the principal values of the D or J tensors 
is much larger than the other two (say J z 3> J x , J y ), spin operators can be 
substituted by two- valued classical variables a p = ±1. In this way, the problem 
reduces to the Ising Hamiltonian 

N N 

H = -J ^2 VpVp+i - VBB^gpdp , (1-17) 
p—i p— l 

in which J and g p may contain information about non-collinearity. In Fig. 1.6 
a sketch of a two-fold, non-collinear Ising chain is shown. Assuming that both 
the g and J tensors have only one non-zero component along their princi- 
pal axes, the corresponding parameters in Hamiltonian (1.17) are given by 

10 Actually, the choice of the unit cell is not unique: one might integrate over any 
pair of internal degrees of freedom dJ?2r- This turns necessary in order to compute 
microscopic averages of individual spin components. 
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J = cos(29)J z and g p — cos(6)g z if B is parallel to the chain axis (c axis in 
Fig. 1.6) while g p = (— l) p sm(8)g z if £? is perpendicular to the chain 11 . The 

Z2r 




slow relaxation fast relaxation 



Fig. 1.6. Sketch of a two-fold non-collinear Ising chain. The local anisotropy axes 
have been chosen coplanar for simplicity and form an angle 6 with the chain axis c. 
For 6 < 7r/4 and J > 0, slow relaxation is expected only when B is applied parallel 
to the chain axis. 

framework in which static properties of SCMs can be modelled seems to be 
well-defined. However, it should not be forgotten that the genuine ID static 
behaviour can be accessed only above a certain temperature Tb, dependent 
on the specific experiment, below which slow dynamics starts playing a ma- 
jor role. Moreover, finite-size effects or 3D interchain interactions may come 
into play at higher temperature than Tb [14]. On the high-temperature side, 
distinctive ID features (short-range order) smear out in the isotropic param- 
agnetic phase. All these phenomena set limitations to the applicability of any 
equilibrium ID model. 

1.5 Glauber model and Single-Chain Magnets 

In this section we will assume the anisotropy at the brick level to be large 
enough that Hamiltonian (1.17) suffices to discuss the important features of 
slow dynamics. A kinetic version of the Ising model was proposed by J. R. 
Glauber in 1963 [44]. As based on stochastic dynamics, this model refers to 
coarse-grained dynamics, typically some orders of magnitude longer than a 
Larmor period. 

Following Glauber, let P(g_, t) be the probability of occurrence of some 
configuration a = <7i, <7 P , un at time t and it; (Jp _j._ (Tp the probability of 
reversing the p-th spin per unit time. The master equation of the problem, 
thus, reads 



11 The reader is addressed to Ref. [38,40] for a more rigorous treatment. 
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^ P fc t) = -^2 w <y p ^-a p P(<?l, , CT p , , CTjv, t) + 

<"« 

+ w -a p ^a p P(vi, ,-<Tp, ,<7N,t) ■ 

p=l 

To model the magnetisation dynamics and a.c. susceptibility it is not necessary 
to solve equation (1.18): one can limit oneself to single-spin averages s p (t) 

s p (t) = ^a p P(a,t) . (1.19) 

It can be shown that spin averages fulfil the differential equation 

d ^ = -2(a pW ^ p ) , (1.20) 

where (...) denotes, again, the time-dependent average performed by means 
of P(a, t). In order that dynamics drives the system towards Boltzmann equi- 
librium, the detailed-balance condition shall hold 

P C q(cri, -Op, ...,<T N ) _ »V-^ ^ 

Peq^l, ...,cr p , ...,a N ) W-„ p -+ ap 

The equilibrium probabilities on the left-hand side of (1.21) are obtained 
reversing the p-th spin while leaving the other (N — 1) unchanged. Their ratio 
can be written as follows 

-Peq(0l,-,-Q"p,--,<Tiv) _ CXp [-K<T p {Op-1 + <7p+l)] CXp (-kpOp) 

P e q(ai,...,a p ,...,a N ) cxp [k<t p {<j p -i + <r p +i)] exp (h k a p ) (122) 
[l — \a p (<7p-i + (Tp+\) tanh(2«)] [1 — a p tanh(/i p )] 
[l + \a p {cf p -i + <J p +i) tanh(2«)] [1 + a p tanh(ft p )] 

with k — J/kT and /i p = ^Bg p /kT. The above relation suggests the follow- 
ing form for the transition probability: 



1 

Wa p ^-a p = 2" 



[1 - <7ptanh(/i p )] , (1.23) 



with 7 = tanh(2«), so that detailed balance (1.21) is automatically fulfilled. 
Equation (1.23) corresponds to Glauber's original choice; other transition 
probabilities fulfilling (1.21) could be chosen [51], but - to our knowledge 
- they have not been considered in the context of SCMs. Note that the pa- 
rameter a entering (1.23) sets the natural time unit of the model. It can be 
interpreted as the attempt frequency of an isolated spin, i.e., the transition 
probability for vanishing exchange coupling, J — 0. Already Suzuki and Kubo 
commented that, in general, a should depend on temperature [52]. We will 
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come back to this important point further on. Combining (1.20) and (1.23), 
a set of differential equations for spin averages is obtained 



1 ds p 
a dt 



1 
2 



7 



(s p -i+s p+ i) + 1 - - (r p _i !P + r P:P+ i) tanh(/i. 



'■)>} 



(1-24) 

where r Pt i = (a p ai). This means that the knowledge of pair-spin correlations 
is needed to solve (1-24). In turn, the knowledge of three-spin correlations is 
needed to obtain r Pi ; and so on. In other words, equation (1.24) is the first one 
of an infinite hierarchy of kinetic equations [44,53]. A judicious truncation of 
this series is, therefore, required in order to get analytic results which could 
easily be compared with experiments. In the following, we will analyse dif- 
ferent decoupling schemes, related to different physically relevant situations. 
Another crucial point concerns the choice of boundary conditions for the sys- 
tem (1.24). A realistic SCM consists of a collection of open arrays of spins 
coupled via the exchange interaction. The length distribution of these arrays 
is determined by the spatial distribution of defects in a sample. In this sense, 
open boundary conditions give a more accurate description of SCM dynamics 
than periodic boundary conditions. However, we start considering periodic 
boundary conditions because calculations are less involved but still provide 
insight into the essential features which are not affected by the presence of 
defects. 



When no external field is applied, h p = 0, the dependence on r p ^ disap- 
pears from (1.24) that then reduces to a linear system of first-order differential 
equations. The corresponding eigenvalue problem involves a circulant matrix 
and is diagonalized by a discrete Fourier transform. A general solution takes 
the form s p = J2 q s q e iqp e~ Xqt , with 

X q = a (1 — 7COsg) (1-25) 

and q — 0, 2n/N, . . . , 2n(N — l)/N set by periodic boundary conditions. The 
initial configuration determines, instead, the Fourier amplitudes. If the sys- 
tem is prepared into a ferromagnetic saturated state with s p = 1 for every p, 
the only nonzero Fourier component corresponds to q = 0, that is So = 1/N. 
Accordingly, the magnetisation is expected to follow a mono-exponential re- 
laxation with a characteristic time scale r = 1/ [a (1 — 7)]. For ferromagnetic 
coupling, J > 0, r diverges exponentially at low temperature like e 4re . Because 
of this divergence, some ferromagnetic ordering may persist over macroscopic 
time scales in the absence of applied field. In this sense, the work of Glauber 
has foreseen what would be observed in SCMs about forty years later. The 
realization of these systems gave the opportunity to generalize the original 
Glauber model to include features of realistic SCMs and specific experiments. 
The response to a tiny a.c. field B = B'e~ lut is modelled by linearising the 
hyperbolic tangent in (1.24). But this does not eliminate the dependence on 
pair-spin correlations. Already Glauber circumvented this problem replacing 
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r p-i,p an d r p,p+i by their equilibrium average, equal to tanh(«;) [44]. Limiting 
himself to equivalent spins (g p — g independent of the site), he predicted that 
the a.c. susceptibility was of the form (1.5), provided that Xeq was taken as the 
static susceptibility of the Ising model and r = 1/ [a (1 — 7)]. As mentioned 
in the previous section, non-collinearity among local anisotropy axes is more 
the rule rather than an exception. This affects Hamiltonian (1.17) through the 
site-dependent Lande factor. The spatial periodicity of g p defines the magnetic 
unit cell. It is worth remarking that the periodicity of g p generally depends on 
the direction along which the magnetic field is applied. The simplest case of a 
two-fold non-collinear Ising chain is sketched in Fig. 1.6. With relatively small 
effort, an analytic formula for the a.c. susceptibility can be derived, which ac- 
counts for non-collinearity or non-equivalence of magnetic centres [54]. For 
lj <C a, a resonant behaviour, i.e. a frequency-dependent peak in x(w), is ex- 
pected only when the field is applied along specific crystallographic directions 
(e.g., the c axis in Fig. 1.6). In particular, those directions are the ones along 
which the ground-state magnetisation is uncompensated. This prediction for 
the dynamic response of non-collinear Ising chains was indeed supported by 
experiment [48,54,55]. 

For several years the truncation schemes summarized above had repre- 
sented the starting point for generalizations of the Glauber model which 
aimed at giving a better account for the characteristics of real SCMs. Then, 
the restriction to zero-field a.c. susceptibility prevented from modelling the 
dependence of relaxation time on static applied field. A breakthrough was 
represented by the work of Coulon and co-workers [56] who actualized the 
local-equilibrium approximation for pair-spin correlations proposed by Huang 
in the seventies [53]. Let us start from refreshing the main ideas of local- 
equilibrium approximation for the case of periodic boundary conditions, as 
treated by Huang. Equivalent magnetic moments, g p — g, coupled ferromag- 
netically will be assumed. With these hypotheses, the single-spin averages are 
independent of the site at thermodynamic equilibrium and read 



m = with Z\ = e- 4K + sinh 2 (^) , (1.20) 



where (. . . ) eq stands for equilibrium average, given by Boltzmann statistics. 
Translation invariance holds also for nearest-neighbour pair-spin correlations, 
of our interest, which are given by 

r-lo a) _ sinh 2 (fe) c~*" [cosh(fe) - 

r (a p+1 a p ) eq — ^— + - [cQsh{h) - A1/2] . (1.27) 

By means of (1.26), sinh(/i) and cosh(ft) appearing in the nearest-neighbour 
correlation can be expressed in terms of m and c~ 4k , which yields 

r = i ^- m ^ (1.28) 

1 + ^m 2 + (1 - to 2 )c 4k 
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In two physically relevant situations translational invariance may be assumed 
for timc-dcpcndcnt spin averages, s p , as well. The first one corresponds to 
having equal initial conditions for all spins: s p (0) = /j,, with — 1 < fi < 1. 
Since g p — g have been assumed, this initial condition is simply realized when 
a magnetic field (possibly zero) has been switched on far in the past [t — > — oo) 
and changed to some different value at time t = 0. The second situation is a 
typical a.c. susceptibility experiment, for which only the stationary response 
to a tiny external drift is relevant. In these two cases, time-dependent spin av- 
erages become independent of the site and the label p can be dropped from the 
variables s p in (1.24). The local-equilibrium approximation consists in assum- 
ing that (1.28), which establishes a closed relation between equilibrium spin 
averages and nearest-neighbour correlations, holds true for time-dependent 
averages as well, namely out of equilibrium. Equation (1.24), thus, simplifies 
as 

1 ds 

- — = - (1 - 7 ) s + (1 - 7 r) tanh(/i) , (1.29) 
a at 

where r is given by (1.28) with m is replaced by s (time-dependent average). 
Within the Glauber model, the local-equilibrium approximation is nothing 
but a trick to truncate the hierarchy of kinetic equations. The resulting equa- 
tions of motion are generally non-linear, the non-linearity arising from r[s]. 
Fortunately enough, equation (1.29) can be solved analytically [53]. More 
importantly, for t — » oo the exact steady-state solution is recovered. For in- 
stance, a mean-field truncation scheme might alternatively be assumed, setting 
r = s 2 , but this would not reproduce the exact steady-state solution. Note 
that local-equilibrium approximation does not require small applied fields. For 
what concerns SCMs, much interest relates to the study of linear departures 
from equilibrium. Following Pini and co-workers [57], let us split the field into 
a static contribution of any intensity (ho = fisBog/kT) plus an oscillating 
field of much smaller intensity B' and with frequency u>: h = ho + ft.'e~ lut . 
As a consequence, s is expected to deviate slightly from its equilibrium value, 
m(T,B ), and (1.29) can be linearised as follows: 

--^ = -(1-7 + 2 7 tanh 2 (fe )) 8s+(l- tanh 2 (/i )) h'e~ iut , (1.30) 

where Ss = s(t) — m and the fact that r[s] w r[m] + (dr/dm) cq 5s with 
(dr/dm) eq = 2tanh(/io) has been used. The stationary behaviour is obtained 

inserting the trial solution Ss — Sse~ lLOt in (1.30), which yields the a.c. sus- 
ceptibility The resulting formula is equivalent to (1.5) and x cq is the suscep- 
tibility that would be obtained by differentiating m in (1.26) with respect to 
B. This matching is a direct consequence of the fact that local-equilibrium 
approximation provides the exact steady-state solution for s. The relaxation 
time is, instead, given by the inverse of the prefactor of Ss in (1.30): 
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The Glauber relaxation time is recovered in the limit ho — >• and, as already 
pointed out, diverges exponentially upon lowering temperature. Note that 
the net effect of a static field is that of removing such a divergence, though 
the dependence of the relaxation time on Bo is much less dramatic than on 
temperature. 

The Glauber model was extended to weakly interacting spin chains by 
Zumer [58]. Similarly to Scalapino [15], he treated the interchain interaction 
as a mean field, limiting his analysis to the critical region around the transi- 
tion to a 3D ordered phase. Equation (1.31) may allow generalizing Zumer's 
results to lower temperatures, away from the critical region. A joint theoretical 
and experimental investigation of this phenomenon would provide important 
information on the critical behaviour of SCMs [14]. A realistic model should, 
however, take into account finite-size effects induced by the presence of defects. 



1.6 Glauber model for finite chains 

Though it may sound somewhat technical, the study of finite-size effects have 
played a central role in theoretical and experimental characterization of SCMs. 
As a first step, open boundary conditions need to be considered, which makes 
the transition probability of extremal spins take the form 

a [1 — t) <7i<72] [1 — o i tanh(/ii)] 

a [1 — rj <jn&n-i\ [1 — &n tanh(Zijv)] , 

with rj = tanh( J/kT), obtained again from the detailed-balance condition. 
The kinetic equations for spin located at boundaries are modified accordingly: 

(si - r\ s 2 ) + (1 - V ri 2 ) tanh(/ii) 

(1.33) 

(s N - t) sat_i) + (1 - r\ rjv-i,iv) tanh(/i jV ) • 

In the absence of external field, the characteristic time scales can be de- 
duced by inserting the trial solution = (A p c lkg + A r c~ %kq ) e~ Xqt into sys- 
tem (1.24) that still holds for bulk spins, with labels 2 < k < N - 1. The 
relation between X q and q remains the same as in (1.25) but the values taken 
by q are different from the case of periodic boundary conditions. Due to the 
loss of translation invariance, both amplitudes A v and A r must be considered. 
A pair of equations for these amplitudes are obtained inserting the trial so- 
lution into (1.33), with \ given by (1.25). For B = 0, this is a homogeneous 
system that only admits the trivial solution A p = A r = unless the determi- 
nant of the coefficients of A p and A r is zero. By requiring this, the following 
implicit equation for the values of q is obtained [59,60]: 



iVcn—t—o-! — 



1 dSi 

a At 
1 ds N 
a dt 
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tan [(N - l)q] = - 



2£ tan q 



(1.34) 



1 - | 2 tan 2 q 



with £ = 77/(7 — v)- The q = solution has to be rejected because it is 



independent of N for every temperature, which is not physical. The remaining 
solutions will be labelled with v, i.e., A„. For ferromagnetic coupling, J > 0, 
the eigenfrequency corresponding to the slowest time scale can be expanded 
for low temperatures to get [61] 



The previous expansion contributed significantly to understanding SCMs. 
From (1.35) one expects the slowest degrees of freedom of the system to equi- 
librate with a relaxation time tn ~ Ne 2K . The fact that the energy barrier at 
the exponent is halved with respect to Glauber's result suggests that, at low 
temperature, relaxation is driven by nucleation of a DW from a boundary. 
At higher temperatures, the Glauber behaviour is recovered. This happens 
when the correlation length becomes significantly smaller than TV and physics 
becomes independent of boundary conditions. Thus, in real systems, the re- 
laxation time is expected to diverge like e 4K at high temperatures - when £ is 
much smaller than the average distance among defects - and like e 2K at low 
temperatures. The experimental observation of such a crossover represented an 
important step in establishing that SCM behaviour could, indeed, be described 
properly in the framework of Glauber dynamics [4,6]. When £ 3> N, the first 
step of relaxation is analogous to the nucleation of a critical droplet to re- 
verse the magnetisation in metallic nanowires or elongated nanoparticles [13]. 
Depending on geometrical characteristics of the sample, non-uniform mag- 
netisation reversal may be favoured with respect to the standard Neel-Brown 
mechanism (uniform rotation). The latter is known to follow an Arrhenius law, 
r ~ e A ^' kT , with an energy barrier proportional to the sample volume. To 
the leading order, the temperature dependence is of the Arrhenius type also 
in the case of non-uniform magnetisation reversal, but A T typically does not 
depend on the sample size. This fact directly originates from the local charac- 
ter of DW excitations which serve as nuclei to initiate magnetisation reversal 
(relaxation), both in metallic nanowires and in SCMs at low temperature. 

After being nucleated at one boundary, a DW may reach the other end 
of the chain with probability <~ 1/N by performing an unbiased random 
walk [27]. This is at the origin of the dependence on N appearing in (1.35) 
and, consequently, in tn- When this is the main channel for relaxation, in real 
SCMs one would expect to observe a decrease of the pre-exponential factor of 
the relaxation time by increasing the number of defects (see (1.3)); the energy 
barrier of the Arrhenius law should, instead, remain constant: A TN = 2J. 
This trend was qualitatively confirmed in Co(hfac) 2 (NITPhOMc) compounds 
in which part of the Co 2+ ions were substituted, in different amounts, by 
non-magnetic Zn 2+ atoms [10,11]. The fact that the pre-exponential factor 




(1.35) 
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increases with the system size is typical of a sizeable time elapsed during DW 
propagation in the relaxation process. In passing, we note that the opposite 
trend, i.e., a decrease of the pre-exponential factor of relaxation time with in- 
creasing the system size, was predicted for nanowires in which magnetisation 
reversal is forced to initiate from the bulk (e.g., in toroidal samples or with 
enhanced anisotropy at the ends) [13]. In that case, the probability to nucleate 
a soliton-antisoliton pair increases with N and the reversal rate consequently. 

The local-equilibrium approximation may also be used to decouple the hi- 
erarchy of Glauber equations when a finite field is applied to an open chain. 
It is convenient to linearise directly (1.24) with respect to 5s p = s p — m p . 
Note that the equilibrium values m p are now site-dependent due to the lack 
of translation invariance. The kinetic equations for Ssk contain the variation 
of nearest-neighbour correlation functions £r p _i iP and 5r p ^ p+ \. For a chain of 
N equivalent spins, with g p = g, Matsubara and co-workers provided a set of 
analytic relations to express equilibrium correlations (a p a p+ i) cq as functions 
of single-spin averages of open chains of different length [62] . If one assumes 
that such relations still hold true out of equilibrium, pair-spin variations can 
be written in terms single-spin averages: 5r p p+1 = A^ p Ss p + B^ p 8s p+ i with 
A NtP and B Np depending only on equilibrium quantities (the reader is ad- 
dressed to Ref. [56] for details). With the same convention introduced in (1.30) 
the response to an a.c. field B'e" 4 "* superimposed to a static field B is de- 
scribed by a system of linear equations of the form: 

— = -MS + a (1 - tanh 2 (fr )) tie'*"* & , (1.36) 

where S = (8si, . . . , Ss N ) T ; the matrix M and the vector \P arc only func- 
tions of equilibrium averages, model parameters, temperature and static field 
(explicit expressions can be found in Ref. [57]). Let (f) u and \ v be the eigen- 
vectors and eigenvalues of M, namely M<j) v = A„0„. The stationary solution 
of (1.36) then reads 

S = a (1 - tanh 2 ^)) h'e~ iut V ^ (1.37) 

\ v 1 — iujT u 

with t„ = 1/A„. The dynamic susceptibility is given by x(w)=<7/iBe ^2 p Ss p /B', 
where Ss p are the components of the £ vector in (1.37). With respect to the 
case with periodic boundary conditions, the choice of a site-independent Lande 
factor does not yield an a.c. response dependent on a single relaxation time. 
The relative weight of different contributions labelled by v shall depend on 
temperature and on the static field Bq. In practice, the matrix M can be di- 
agonalizcd numerically for any values of B and T. The size of this matrix, N 
by N, is set by the number of spins in the chain. Realistic values of N fall in 
the range 10— 10 4 , meaning that x(w) can easily be computed with standard 
diagonalisation routines. Among other things, this allows checking whether a 
unique relaxation time is dominating the summation (1.37) and thus x(w). 



1 Single-Chain Magnets 25 



When the distribution of defects in a SCM compound is not peaked, an aver- 
age over all the possible lengths may be required to compare the theoretical 
susceptibility with experiments [10,11,57]. Analogously, for a comparison with 
experiments on powder samples an average over all the possible orientations 
of the applied field with respect to the easy axis might be needed [56]. Due 
to space limitations, we prefer not to enter the details of those averaging 
procedures but, rather, address to the existing literature. 

The divergence of relaxation time upon lowering temperature can be inter- 
preted as critical slowing down. The ID Ising model displays a magnetic phase 
transition at zero temperature, meaning that the critical point is located at 
the origin of the (T, Bo) plane, that is T = and Bq = 0. Since the divergence 
of the correlation length is hampered by defects, it is more appropriate to in- 
vestigate the critical behaviour of SCMs with finite-size scaling. For B = 0, 
Luscombc et al. noted that the ratio between the relaxation time of a finite 
chain, tn, and that of the infinite chain, r originally obtained by Glauber, is 
a universal function of x = N/£, when both N, £ »1: 

T = f{x) = /m\ 2 ' (L38) 

where u(x) is the smallest root of the transcendental equation wtan(w/2) = 
x [60]. By definition, f{x) tends to one for x » 1; while for x <C 1 it is 
f(x) ~ x/2. The reader may easily verify this limit by using formulae t ~ 
c 4K /2a, t n ~ Ne 2K /2a and £ ~ e 2K /2, which hold for TV, £ >1 (see (1.35) 
and (1.25)). More recently, Glauber dynamics of the open chain in presence 
of realistic fields was studied by Coulon and co-workers [56] who found 

T N (B Q = 0) _ 1 , 2 2 



tat (B ^0) 



= l + a 2 /io; (1-39) 



remarkably, the constant on the right-hand side is given by a = 2£/(x^/2/3), 
f(x) being the scaling function defined in (1.38). For x ~> 1, the limit a = 
2£ is recovered by expanding the hyperbolic tangent in (1.31) (remember 
that f(x) — > 1 in this limit). In the opposite limit, one has a = y // 2/3N 
consistently with the work of Schwarz developed in the context of helix-coil 
transition of polypeptides [63]. The quadratic dependence on Bo of the ratio 
of relaxation times in the vicinity of the critical point stated by (1.39) was 
confirmed by experiment: first, in SCMs made up of repeating trinuclear units, 
Mn 3 +-Fe 3+ -Mn 3 + and Mn 3 +-Ni 2 +-Mn 3 + [56], later in Co(hfac) 2 (NITPhOMe) 
compounds [57] . As pointed out in Ref . [56] , the quadratic dependence on B is 
also expected for SMMs. In fact, when repeating units in a spin chain consist of 
SMM-like centres an additional dependence on temperature and on B enters 
the Glauber model through the attempt frequency a. Thus, information about 
the ID universality class is somehow contained in the scaling function f(x) 
rather than in the quadratic take-off of tn (Bo ^ 0) as a function of the applied 
field. 
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In summary, the Glauber model prescribes that precise relations among 
characteristic energy scales shall hold for a text-book SCM. Recalling (1.3) 
and (1.4), the barrier controlling the divergence of the correlation length can 
be directly deduced from static susceptibility measurements at high enough 
temperature. The last condition is required in order for £ to be smaller than 
the distance among defects. When this ceases to hold true, a saturation of 
the product is observed, at low T. According to the Ising model A^ = 
2J, which implies that the energy barrier of the relaxation time is expected 
to be A T = 2A^ = 4J at high temperature and A Tn = A^ = 2J at low 
temperature. Moreover, the crossover between the thermodynamic limit (£ <C 
N) and the finite-size regime (£ 3> N) should be described by (1.38) and (1.39) 
in the presence of a static applied field. SCMs represent a class of model 
systems in which most of these predictions were confirmed. Often, finding a 
quantitative agreement required ad-hoc generalizations of Glauber's idealized 
picture, without renouncing its basic concepts. Some of those generalizations 
will be discussed in the next section. 



1.7 Beyond the Glauber model 

An important generalisation of the Glauber model relates the temperature 
dependence of the parameter a [64]. Introduced in (1.23) for the single-spin 
transition probability, this parameter turns out to be the proportionality coef- 
ficient between the low-temperature expansions of the correlation length and 
the relaxation time: r = 2£ 2 /a. In other words, for the time scales of interest, 
one has that a — 4D S , with D s being the diffusion coefficient for thermally- 
driven DW motion (see (1.2)). Given this equivalence, we will focus on D s 
henceforth. For explaining the experimental results of a SCM made of Mn 3+ - 
Ni 2+ -Mn 3+ repeating units it was proposed that D s ~ e -^A/kT^ wnere /\ A 
was the global effective anisotropy energy of each unit [4]. The relationship 
between the energy barriers of the correlation length and the relaxation time 
was adapted accordingly: A T = 2A^ + A a- The last formula has been vali- 
dated by experiments on a variety of SCMs with sharp DWs. In those cases, 
it was also found that the energy barrier of r has to be modified consistently 
at low temperature, namely A TN = A^ + Ax- 

One minor remark is that the Z\a contribution to the energy barrier of the 
relaxation time is justified only when some single-ion anisotropy is present. 
For instance, we have seen that Co 2+ in distorted octahedral environment can 
be assumed to behave as an effective spin one-half at low temperature. This 
assumption is not consistent with a finite A\ for SCMs based on Co 2+ . More 
importantly, the picture appears more blurred for broad DWs. Let us refer 
again to Hamiltonian (1.1). As mentioned in Sect. 1.2, for J > D the corre- 
lation length in units of DW width is a universal function of the temperature 
expressed in units of DW energy: £/w = A (S^/kT). It has been known since 
the eighties that a spin wave can propagate across a broad DW acquiring a 
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phase shift 12 [65]. In order to conserve the total magnetisation at short time 
scales, the DW is left displaced after this scattering event [66]. When many of 
such events occur incoherently and involve thermalised spin waves, the result- 
ing DW motion may be assimilated with that of a Brownian particle. Indeed, 
still in the eighties, it was shown that D s scales like the square of the ratio 
kT/Sdw in the absence of damping [67] and linearly when some damping term 
is included [68] . Dimensional analysis suggests to complete the latter result as 



•• 2 kT 



D s a— — , (1.40) 

Td Odw 

with Td being a characteristic time scale of the problem, associated with short- 
time dynamics. Equation (1.40) is consistent with recent numerical results re- 
ported in Ref. [30] . In the same paper, the activated behaviour of D s expected 
for sharp DWs was recovered as well. A qualitative argument for the different 
temperature dependence of D s expected for sharp and broad DWs can be 
given starting from zero-temperature dynamics. In the continuum formalism 
one finds that a field of any intensity applied along the easy axis is able to 
move a broad DW [24,65,70,71]. In the opposite limit, it was shown that 
a finite threshold field is needed to let a sharp DW propagate [26]. In this 
case, translating a DW requires local modifications of the spin profile, which 
create an effective Peierls potential. This potential is periodic with respect 
to the position of the DW centre and the difference between its minima and 
maxima decreases exponentially with increasing the DW width [28,72], till 
it vanishes in the continuum limit. It seems, therefore, plausible to expect a 
thermally-activated diffusion coefficient only for sharp DWs. 

While for sharp DWs the relaxation time depends on J, D and T indepen- 
dently, our present understanding of SCMs suggests that r should depend only 
on the ratio £& w /kT for broad DWs. This can be readily deduced by relating 
r to the correlation length £ = wA (S^/kT) by means of the random- walk ar- 
gument and (1-40) (remember that this argument holds only for £ <C N) [30]. 

The standard theoretical framework to deal with magnetisation dynam- 
ics is the Landau-Lifshitz-Gilbert (LLG) equation. In that context, one ex- 
pects Td introduced in (1.40) to be of the order of the dumping time: 
r d ~ (1 + a G )/(aG7o#A) — ^/(2-Dac), where qq <C 1 is the Gilbert damp- 
ing [73], 70 is the gyromagnetic factor and Ha the anisotropy field. For values 
of D that are realistic for SCMs, h/D falls in the picosecond range while 
the damping constant is typically aG=10 _1 -10 -4 . As slow dynamics is usu- 
ally probed in SCMs at time scales larger than milliseconds, clearly it per- 
tains to long-time behaviour in the language of LLG equation. Moreover, 
since physics of SCMs is dictated by thermal fluctuations, a stochastic noise 
should be included in numerical simulations [74-76] . In spite of the enormous 
improvements experienced in computational capabilities [77], performing a 



12 Recently, magnonic applications of DWs which exploit such a phase shift has been 
proposed in the context of metallic nanowires [69]. 
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stochastic-dynamic simulation which covers a time window of several orders 
of magnitudes still remains prohibitive. In this sense, the brute-force approach 
to SCMs dynamics does not seem promising for the next future. 

With respect to the sharp-wall case, there is no analogous of the Glauber's 
formalism for SCMs with broad DWs. Experimental realizations basically con- 
sist of ferrimagnetic chains alternating Mn 3+ with an organic radical [43] . A 
reasonable model is the one which produces the kernel (1.15), where Mn spins 
are treated as classical vectors. In the experimentally accessible region A^ can 
be much smaller than £dw _ up to about half of it - due to spin-waves renor- 
malisation [9]; while for sharp DWs one has A^ = £dw at any temperature 
kT < J [30]. This fact needs to be taken into account in the experimental 
characterization of SCMs with broad DWs (see Fig. 1.2). For what concerns 
the barrier of the relaxation time, the available experimental results yield A T 
about 10-20% times larger than £dw 13 [9,43]: much smaller than twice the 
DW energy at T=0 as predicted by Glauber for the Ising model. Making a 
definitive statement about the origin of energy scales involved in dynamics is 
not possible yet. Defects probably affect the nucleation and diffusion of broad 
DWs differently with respect to the Ising limit. In metallic nanowires, for in- 
stance, defects act as pinning centres for DWs or vortices. In SCMs a similar 
phenomenon may induce a reduction of DW mobility, namely D s . Another 
possibility is that DWs may preferentially be nucleated at defects because it 
is energetically favourable 14 . Only a thorough characterisation of SCMs with 
broad DWs in which the concentration of defects may be controlled could al- 
low answering those questions. At the same time, such a study would provide 
important information about the joint effect of defects and thermal fluctu- 
ations. This would also be relevant for DW dynamics in metallic nanowires 
that are typically described by the very same classical Heisenberg Hamilto- 
nian (1.1) [13,28,71,78-81]. 

1.8 Conclusion and perspectives 

The title of the review we wrote about five years ago was "Single-chain mag- 
nets: where to from here?" [37]. The idea was that of reviewing critically 
what had been done in the synthetic, experimental and theoretical fields. The 
analysis indicated that the hunt for high-temperature blocking magnets was 
going to continue. This has been confirmed but with the explosion of the 
interest for Lanthanides with the challenging difficulties associated with the 
large unquenched orbital moment [82]. Much more work shall be done espe- 
cially in theory. Another field which is developing fast is that of ab initio, 



For these mixed chains £dw = 2VJD, with a factor 2 of difference with respect 
to Hamiltonian (1.1) [9]. 

Accommodating the DW centre onto a defect reduces the anisotropy energy and 
£dw consequently. 
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DFT calculations which are rapidly complementing/substituting Ligand-ficld 
approaches [47,83]. Far-from-equilibrium dynamics and aging [18,84] as well 
as the interplay between SCM behaviour and quantum effects [19,85] call for a 
more systematic investigation. The comparison of the properties of molecular 
nanomagncts with elongated magnetic nanoparticles and magnetic nanowircs 
has been stated a few times throughout the chapter [13,30]. We feel that 
SCMs can provide good insight into the finite-temperature behaviour of such 
nanosystems. Finally, as a matter of facts, molecular systems have already en- 
tered the domains of spintronics [86] and quantum computing [87] . In future, 
besides their traditional role as model systems, SCMs can possibly find their 
place in those applicative research contexts. 
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